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ABSTRACT 

We report on multi-frequency, wideband radio observations of the Galactic Center magnetar 
(SGR 1745—2900) with the Green Bank Telescope for ^100 days immediately following its initial 
X-ray outburst in April 2013. We made multiple simultaneous observations at 1.5, 2.0, and 8.9 GHz, 
allowing us to examine the magnetar’s flux evolution, radio spectrum, and interstellar medium pa¬ 
rameters (such as the dispersion measure (DM), the scattering timescale and its index). During two 
epochs, we have simultaneous observations from the Chandra X-ray Observatory, which permitted 
the absolute alignment of the radio and X-ray profiles. As with the two other radio magnetars with 
published alignments, the radio profile lies within the broad peak of the X-ray profile, preceding the 
X-ray profile maximum by ^ 0.2 rotations. We also find that the radio spectral index 7 is significantly 
negative between ^2 and 9 GHz; during the final ^30 days of our observations 7 ^ —1.4, which 
is typical of canonical pulsars. The radio flux has not decreased during this outburst, whereas the 
long-term trends in the other radio magnet ars show concomitant fading of the radio and X-ray fluxes. 
Finally, our wideband measurements of the DMs taken in adjacent frequency bands in tandem are 
stochastically inconsistent with one another. Based on recent theoretical predictions, we consider the 
possibility that the dispersion measure is frequency-dependent. Despite having several properties in 
common with the other radio magnetars, such as Tx,qui/Arot ^ 1 ? an increase in the radio flux during 
the X-ray flux decay has not been observed thus far in other systems. 

Subject headings: Galaxy: center — pulsars: individual (PSR J1745—2900, SGR 1745—2900) — stars: 
magnetars 


1 . INTRODUGTION 

Magnetars are exotica among the exotic: whereas other 
pulsars are sustained by their stored angular momentum, 
the primary energy source that powers this special class 
of ob jects is likely the neutr on star’s immense magnetic 
field (jMereghetti et al.ll2015f ) . The field strengths take on 
the highest values ever inferred, typically > 10 ^^ G and 
even up to ^ 10 ^^ G. According to the McGill Online 
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Magnetar GataloeF^ (|01ansen fc Kaspill2Q14D . there are 
28 known magnetars, of which only four have displayed 
pulsed radio emission. 

SGR 1745-2900 (J1745-2900, hereafter) is the most 
recent addition to the small collection of magnetars with 
observed pulsed radio emission (the “radio magnetars”, 
to which we will refer by their PSR names: J1809—1943 
(XTE 1810- 197), J1550-5418 ( I E 1547 . 0-5408), fc 
J1622—4950 (|Camilo et al.l 120061 l2007al: iLevin et al.l 
Imoh ). On 25 April 2013, one day after the XRT 
aboard the Swift satellite detecte d flaring activity coin - 
cident with the Galactic Genter (jPegenaar et al.ll2013[ ). 
a short X-ray burst was observed by Swift/BAT show¬ 
ing characteristics similar to those usually obse rved from 
soft gamma-ray repeaters (jKennea et al.ll^l3c[ ). Shortly 
thereafter, observations from the NuSTAR satellite iden¬ 
tified the source as a magnetar with a Ps = 3.76 s spin 
period, and its radio pulsations we re subsequently se en 
by the Effelsberg lOQ -m Telescope (|Mori et al.ll2013allbl : 
lEatongh et al.l[^13af ). J1745—2900 was soon physically 
associated with the Galactic Genter, located only ^2.5” 
away from Sagittarius A* (Sgr A*) with a neutral hydro¬ 
gen column density and dispersion measure (DM) consis¬ 
tent with b eing within ^2 pc of the Milky Way’s central 
black hole (jEatongh et al.l[2013bl: iRna et al']l2013f ). 

Early determinations of its spin-down Pg put 
J1745—2900 squarely within the magnetar population, 
having an inferred magnetic field strength at the equator 

http://www.physics.mcgill.ca/-pulsar/magnetar/main. 
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Bs ^ 3.2 X 10^^ G ^ 1.6 X 10^"^ G, a characteristic 

age Tc ^ Ps/{2Ps) ^ 9 kyr, and a spin-down luminos¬ 
ity of ^ = Lrot = 3.95 X 10^^ erg s~^{P~^Ps) ^ 4.9 x 
IqSS g-i (jPea et al.ll2Q13f ). However, its estimated 

J uiescent X-ray lumino sity of Lx,qui < 10^^ erg s“^ 
Coti Zelati et al.l[2Q15[ ) may place J1745—2900 on the 
side of I/x,qui/^rot < 1? opposite the “classic magne- 
tars” but alongside the other three radio magnetars, 
high-H pulsars, and radio pulsars with X-ray emission 
(IRea et, a,lJl2iTI^ . 

Given the unique environment in which J1745—2900 
resides, the detection of its radio pulses is some¬ 
what surprising. Indeed, numerous surveys of the 
Galactic Genter region covering ^1-20 GHz have 


failed to find a pulsar within the central i 

Darsec 

(most recentlv. iJohnston et al.ll2006 

: iDeneva et al.ll2009l: 

Macouart et al.l I2OIOI: iBates et al.l 

201 It ISiemion 

et al.l 

20131). The discovery of this single magnetar has 


l ed to a windfall of implications fo r future discoveries 
(iGhennamangalam fc Lorimerl I20l4 iDexter fc O’Lear^ 
|2014 iMacQuart fc Kanekarll2015[) . Because of its prox¬ 


imity to the Galactic Genter, J1745—2900 has the largest 
DM (1778 cm“^ pc) and rotation measure (—6.696 x 
10^ ra d m“^) of any known pulsar (|Eatongh et al.l 
l2013b[ ). The predicted value for the scattering 
timescale at 1 GHz, based on empirical relationships 
given its DM, is ^ 1000 s (|Krishnaknmar et al.l 120151: 
iLewandowski et aH l2015aD . meaning that J1745—2900 
would be undetectable at frequencies less than ^5 GHz. 
The situation is exacerbated by the presence of an 
additional scattering screen in the Galactic Genter 
(|Gordes fc Laziol I1997D . Normally, the prospect of de¬ 
tecting distant radio pulsars above several G Hz is bleak, 
since their average spectral index is ^ -1.4 (|Bates et al l 
I2013D . However, because the other radio magnet ars 
have flat/inverted spectra, one might expect to detect 
J1745—2900’s unscattered pulse profile at high frequen¬ 
cies. In the analyses that follow, we will reiterate the 
finding that J1745—2900 has a s ignificantly smaller scat¬ 
tering timescale than predicted (jSpitler et al.ll20l4 ). and 
will show that J1745—2900 was much brighter at lower 
frequencies, having a very negative spectral index some 
100 days after the onset of its outburst , even though 
more recent observations by iTorne et~an (j2015[ ) showed 
the spectral index has since flattened. 

In this paper, we analyze multi-frequency radio 
data over the first ^100 days after J1745—2900’s 
discovery, during which time there were two additional 
Swift / BAT-detected bursts on 7 June 2013 and 5 August 
2013 (jKennea et al.l [20T3a]lbl ). For two of our epochs, 
which bracket the third burst by ^1 week on either side, 
we have simultaneous Chandra observations. These 
observations allow us to find the absolute alignment of 
the radio and X-ray profiles, and to look for correlated 
events. We comment on the spin evolution and timing, 
and examine the profile stability, the radio flux evolu¬ 
tion, and the radio spectrum. Finally, we make global 
models of the profile evolution across the low frequency 
bands in order to examine the temporal and frequency 
dependencies of the scattering timescale and dispersion 
measure. We then discuss characteristics of this source 
in comparison with other radio-loud magnetars. 


2. OBSERVATIONS 
2.1. Radio 

We made early detections of J1745—2900 during four¬ 
teen observing epochs with the 100-m Robert G. Byrd 
Green Bank Telescope (GBT) in three different frequency 
bands with various overlap: 1.1-1.9 GHz (5 epochs), 1.6- 
2.4 GHz (7 epochs), and 8.5-9.3 GHz (11 epochs) (PI: 
A. Possenti). Because each observation covers a large 
bandwidth, we refer to each set of data based on the 
IEEE radio band for which each of the receiver systems is 
named (“L-band”, “S-band”, or “X-band”, respectively), 
instead of referring to specific (central) frequencies. Ta¬ 
ble [1] contains details of the observations. In all cases, we 
observed using the Green Bank Ultimate Pnlsa r Process¬ 
ing Instrument fGUPP f^JDnPlain et al.ll2QQ8[) in “inco¬ 
herent search mode”, recording dual-polarization time- 
series data in 2048 frequency channels with a temporal 
resolution of 0.65536 ms. 

Each epoch’s data were folded with the pul¬ 
sar software library DSPSrF^ using a nomi¬ 
nal ephemeris with a constant spin frequency 
(see ^3.3[) and the Chandra-deteimined position 

g.T^nnnn = Sj2000.0 = 29°00'29'/84 

(|Rea et al.ll2013[) . The data were initially folded into 1 
min subintegrations, with 2048 profile phase bins across 
128 frequency channels. We adopted the published 
dispersion measure value of 177 8 cm~^ pc for averagin g 
frequency channels together (jEatough et al.l l2013bf ). 
Persistent, narrow-band radio frequency interference 
(RFI) was excised automatically; any remaining sig¬ 
nificantly corrupted channels or subintegrations were 
removed from the data by hand. 

Galibration scans were taken for each observation us¬ 
ing the local noise diode, pulsed at 25 Hz while on 
source. We recorded on- and off-source scans of a stan¬ 
dard flux calibrator (QSO B1442+101) in each frequency 
band only during the final epoch (MJD 56516). We 
have used this one set of flux calibration scans to cal¬ 
ibrate the whole data set. Standard pr ograms from 
the PSRCHIVEF^ pulsar sof tware library (|Hotan et al.l 
I2QQ4 Ivan Straten et al.ll2Q12[ ) were used to calibrate the 
absolute flux density scale of the noise diode, which is 
then used to determine the magnetar’s flux densit\F1. 

The combination of the large amount of observed scat¬ 
tering f li3.5.ip . the pulsar’s spectrum f 113.4.31 Figure [7|), 
receiver roll-off, and the presence of gain variations (see 
below) rendered significant portions of the ends of L- 
band useless. Namely, there was no pulsed signal in the 
lower 300 MHz portion of L-band, which we masked from 
further analysis, along with the top 50 MHz (which is 
part of the overlap with S-band). In combination with 

14 yyy ^ Safe. iirao.edu/wiki/bin/view/CICADA/NGNPP 
http://dspsr.sourceforge.net/ 
http://psrchive.sourceforge.net/ 

The PSRCHIVE calibration process produced unphysical re¬ 
sults for the earliest S-band detection (MJD 56424); we have cal¬ 
ibrated it by using an approximation based on the measured S- 
band system e quivalent flux densi t y and the radiometer equation 
(cf. §7.3.2 of ILorimer &: Krameii ((20041 b. The result is reason¬ 
able, given that the next S-band observation five days later has a 
comparable fiux density (see Figure [G]). 
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the narrow-band RFI, this left less than ~400 MHz of 
clean, usable bandwidth. Similarly, at S-band we had to 
remove the lower ^100 MHz and the upper ^25 MHz, 
and in total ^625 MHz of usable band remainecQ. Only 
3% of the data was clipped from either end of X-band, 
with a total of 10% removed. We took these seemingly 
draconian measures to offset the original data quality and 
to ensure that the time- and frequency-averaged profiles 
were of reasonably high quality (e.g., see Figured]). This 
was enabled by the source’s relatively large flux density. 

The data quality situation at X-band was s till m ore 
complicated. As also noted by iLvnch et al.l (j2Q14f ) in 
their investigation of this magnetar, large gain variations 
on timescales from a fraction of a pulse period to several 
seconds (visible in the time-series data) are prevalent in 
X-band at the GET, when pointed at the Galactic Cen¬ 
ter. The variations did not (necessarily) integrate away 
over hour-long observations and are representative of a 
stochastic red-noise process. We attribut e these vari¬ 
ation s to changes in atmospheric opacity (jLvnch et al.l 
l2Q14f ) and/or small pointing errors, noting a strong res¬ 
onance in the GET X-band pointing very near 0.3 hS 
The gain variations would be manifested by the relatively 
small beam of X-band (^1.4', compared to ^6' and ^9' 
for S- and L-band) oscillating over the crowded, bright 
Galactic Center (the central parsec extends ^0.4', and 
the separation o f J1745—2900 from Sgr A* is only ^0.04' 
(|Rna et al.ll2013f )). Additionally, it is likely that the base¬ 
line variations are much less prominent at low frequencies 
because they act as “zero-DM” signals that get smeared 
out w hen the pulsar’s signal is dedispersed. iLvnch et al.l 
(|2014f ) also state that the effect may be a function of el¬ 
evation, which fits with our pointing-resonance hypothe¬ 
sis, since the influence of variable elements like the wind 
will be a function of elevation. The persistence and vari¬ 
ability of these variations can be seen in Figure [2l 

The analyses that follow utilized these folded profiles 
in a variety of reduced forms. Unless otherwise noted, 
the reduced radio data have 2048 profile bins (^7.2 ms 
per bin), 32 frequency channels (25 MHz per channel), 
and 5 min subintegrations; in this work, we only consider 
the total intensity profiles. 

2.2. X-ray 

During two of our radio epochs, MJD 56500 and 
MJD 56516, we obtained simultaneous observations of 
J1745—2900 with the Chandra X-ray Observatory (Obs. 
IDs 15041 & 15042; PI: D. Haggard). Table [2] contains 
details of the X-ray ob servations (for further details see 
iCoti Zelati et al.ir2015f ). The field of the first observation 
is shown in Figure [2 the second observation was essen¬ 
tially the same. In each observation, J1745—2900 was po¬ 
sitioned on the back-illuminated chi p S3 of the Advanced 
GGD Imaging Spectrometer (ACIS. lGarmire et al.ll200^ 
instrument. The data were reprocessed with the Chan- 

In two epochs, however, instrument problems left only half of 
S-band viable. See Table [T] 

Even though the average pointing errors at X-band are only 
on the order of several arcseconds at mid-elevations and mild wind 
conditions, the power spectrum in elevation offset shows reso¬ 
nances overlapping with the magnetar’s spin frequency (0.27... Hz). 
See http :// www.gb.nrao.edu/~rmaddale/GBT/Commissioning/ 
Pointing_Gregorian_HighFreq/PntStabilityXBand.pdf for 

details. 


Table 1 

Summary of GBT Observations 


UTC 

MJD 

Bands 

Approx. Length 

Epoch 


Observed 

[min] 

2013-05-04 

56416 

X 

20 

2013-05-12 

56424 

s*,x 

122,200 

2013-05-13 

56425 

X 

60 

2013-05-14 

56426 

X 

49 

2013-05-17 

56429 

s,x 

70,53 

2013-05-23 

56435 

X 

50 

2013-05-30 

56442 

X 

58 

2013-06-21 

56464 

X 

54 

2013-07-14 

56487 

X 

71 

2013-07-15 

56488 

L,S 

120,132 

2013-07-27 

56500 

L,S,X 

186,108,68 

2013-07-28 

56501 

L,S* 

133,117 

2013-08-03 

56507 

L,S 

112,75 

2013-08-12 

56516 

L,S,X 

120,60,56 


Note. — The listed dates and MJDs for the epochs 
are representative of the majority of the epoch, not the 
start time; observations on the same day were taken in 
tandem. The two boldfaced epochs are those for which 
we have simultaneous observations with Chandra. The 
lower half (400 MHz) of the two S-band observations 
with an asterisk were corrupted and unusable. The 
horizontal lines separate the epochs during which the 
three observed types of X-band profile are seen (see 
El and Figure!^. 
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Figure 1. Examples of L- and S- band profiles averaged over all 
epochs. The profiles are shown with 1024 phase bins for clarity. 
Th ese d ata are aligned via a wideband portrait model, as described 
in 3Ml In general, the un-averaged profiles were also of good 
quality, with only minor systematics in the baseline. The total 
bandwidth covered across these two bands is about 1 GHz, from 
~1.4 to 2.4 GHz; 25 MHz of data were averaged for each of these 
profiles, with their center frequencies shown. The profiles were 
very well described by a single scattered Gaussian component, and 
so we do not over-plot the wideband model. The vertical dotted 
lines show example s of o n-pulse regions used for the flux density 
measurements. See E31 for details. 
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Figure 2. Examples of time- and frequency-averaged X-band pro¬ 
files. The profiles are shown with 1024 phase bins for clarity. The 
baseline variations were removed on a profile-to-profile basis by fit¬ 
ting a high degree polynomial (red dashed lines) to the off-pulse 
region (outside the do tted lines) in order to make measurements of 
the flux density (see El- The on-pulse phase window varied in 
size between about 6 and 8%. T he p rofile evolved monotonically 
from one “type” to the next (see ^3.21 and Table [T]). 


dra Interactive Analy sis of Observations so ftware pack¬ 
age (ciAO, version 4.6. lFrnsciQne et aI]l2QQ6f ) and the cal¬ 
ibration files in the caldb release 4.5.9. 

In both observations, J1745—2900 was bright enough 
to cause pile-up in the ACIS detector. A “pile-up map” 
created with the Cl AO tool pileupjnap confirmed that 
mild pile-up was present. Exclusion of data near the 
center of the point-spread function (PSF) from the anal¬ 
ysis would have resulted in the loss of too many photons 
(63% of the source counts were in the two central pixels). 
Moreover, the external part of the PSF contained a sub¬ 
stantial number of counts from Sgr A*. We thus decided 
to proceed as follows. 

We extracted the source counts from a circular region 
centred on J1745—2900 with a 1.5" radius (see Figure [3j); 
this region includes the piled-up events. This area cov¬ 
ers ^85% of the Chandra PSF (encircled energy frac¬ 
tion) at 4.5 keV. A larger radius of 2-2.5" would let in 
more counts from Sgr A* and would only marginally in¬ 
crease the encircled energy fraction. Because of the com¬ 
plex environment, the background spectrum needed to 
be extracted close to the source. We used a thin an¬ 
nulus (with radii of 2" and 4"), excluding a bright area 
associated with Sgr A*. The spectra, the ancillary re¬ 
sponse files and the spectral redistribu tion matrices wer e 
created using specextract. Following iRea et al.l (|2Q13f ). 
we adopt a pure blackbody for the spectral shap e. We 
correc ted the spectra using the pile-up model by iDavisI 
(|2QQ1[) , as im plemented in the mo deling and fitting pack¬ 
age SHERPA ((Freeman et al.|[2QQTI) . The pile-up fraction. 


Table 2 

Summary of Simultaneous Chandra Observations 


Obs. ID 

Radio epoch Exposure time 

Net source 

RMS pulsed 


[MJD] 

[ks] 

counts [10^] 

fraction [%] 

15041 

56500 

45.4 

15.7 

28.8 ± 1.5 

15042 

56516 

45.7 

14.4 

28.9 ± 1.8 

Note. 

^ The 

Icr uncertainties for the 

RMS pulsed 


fractions were deter mined from Monte Carlo simulations (cf. 
iGotthelf et al.l (|1999lP . By another measure, the pulsed fractions 
— defined as the difference between the profile maximum and 
minimum divided by their sum — are ~48%. The folded profiles 
are shown in Figure 



Figure 3. Chandra field of J1745—2900 for observation 15041. 
1 ACIS pixel = 0.492''. The source counts were taken from the 
central-most encircled region (red circle). Background counts were 
extracted from the annulus between the outer two (yellow) circles, 
excluding th e are a marked as “Sgr A*”. We account for pile-up as 
described in G21 


estimated by fitting the jdpileup model, is 3.7% for the 
first observation, and 4.1% for the second. We did not 
attempt any correction of the light curves; the pile-up 
fraction is modest and, in general, pile-up affects sp ectra 
more than it does light curves and pulse profilefd. The 
spectral model fits were acceptable only when the pile- 
up model component was included. A summary of the 
spectral fits is given in Table [3l 

3. RESULTS 

3.1. Transient Events 

J1745—2900 is kno w n to show nar r ow ind i vidual pulses 
(iSoitler et al.l l20l4 iBower et al.l l20l4 iLvnch et al.l 
1201 41). similar to the radio magnetar J1622—4950 
(|Levin et al.l l2012f) . We performed a cursory analysis 
of J1745—2900’s individual pulses in our X-band data, 
seeking only to find anomalous burst-like events in the 
radio data that might be coincident or correlated with 
X-ray features or flares. For this, we took two ap¬ 
proaches. In the first case, we folded the raw data 
into single-rotation integrations, approximately main¬ 
taining the original temporal resolution, averaging over 

This is true unless the pulse profiles are strongly dependent 
on energy, w hich is not the case for J1745—2900, though we refer 
the reader to lCoti Zelati etldl (|2015l i for further details. 
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Table 3 

Chandra Spectral Results 


Obs. ID 


r 

A^h 

kT^ 


Observed fiux*^ 

Luminosity‘s 

Xred (dof) 




[10^^ cm“^] 

[keV] 

[km] 

[10“^^ erg cm“^ s“b 

[10^^ erg s“b 

15041 

0-50io.05 

99.8% 

1.26 ± 0.03 

0.82 ± 0.01 


8.9± 1.2 

o 9 + 0.4 
'^•^-0.3 

1.00 (288) 

15042 

0-48tS:?l 

97.1% 

1.23 ± 0.03 

0.83 ± 0.01 

2 IR + O-13 
^•-LO_o.16 


2.8 ± 0.4 

1.00 (287) 


Note. — The abundances used in the absorbed blackbody model are th ose of I Anders &: Grevessel (| 19891) and photoelectric 
absorption cross-sections are from lBalucinska-Church &: McCammonl (|1992r) . See lCoti Zelati et a D (l2Q15r) for a complete treat¬ 
ment of these observations in the context of a long-term X-ray monitoring campaign. Parameter uncertainties in the table are 
la. 

^ Parameters of the jdpileup SHERPA pile-up model; /x is the grade -migration p arameter and / is the fraction of the PSF treated 
for pile-up, required to be in the range 85-100%. For details, see lDavil (I2001f) and “The Chandra ABC Guide to Pileup”. 

^ The blackbody temperature and radius are calculated at infinity and assuming D = 8.3 kpc (|Genzel et al.|[2010r) . which is 
assumed throughout this work. 

In the 0.3-8 keV energy range; for the luminosity we again assumed D = 8.3 kpc. 


frequency, and summing the polarizations. These data 
were inspected visually. In the second case, we an¬ 
alyzed the raw data with the PrestcFI pulsar soft¬ 
ware package. Here, we applied an RFI mask to 
the raw data with rfifind. We then made a dedis- 
persec£3, frequency-averaged time-series with prepdata 
for each X-band epoch, and searched for single pulses 
with the boxcar-convolution algorithm implemented in 
single_pulse_search.py. We repeated this process on 
the unmasked raw data. 

Single pulses were detected; indeed, one to several 
pulses are visible by eye during almost every rotation. 
However, we saw no anomalously large single pulses or 
other bursts. The distributions of estimated single pulse 
energies all peak at <1 times the average profile energy 
and were inconsistent with power-law distributions. The 
phases of the single pulse arrival times were consistent 
with occuring within the on-pulse window, and the dis¬ 
tributions of the resolved single pulse widths peaked near 
3-4 samples ^ 2 ms, in agr eement with t he X- band scat¬ 
tering timescales found by iBower et al.l (IMl . 

Similarly, the (unfolded) X-ray light curves during the 
two simultaneous observations, binned from 0.5 to 5000 s, 
were featureless and constant. The probability of 
constancy was high for both observations, regardless of 
the choice of binning (>30% and frequently approach¬ 
ing 100%). Due to the uniformly poor quality of the 
X-band data as previously described, we refrain from fur¬ 
ther analysis or discussion of this aspect of J1745—2900 
and direct the reader to the observations of its X-band 
single pulse s as observed with t he V ery Large Ar r ay an d 
the GET in IBower et al.l (|2014D and lLvnch et al.l (|2014D . 
respectively. 

3.2. Profile Variability 

Figure [2] shows examples of the three general types of 
observed X-band profiles, as well as corresponding exam¬ 
ples of our baseline removal and on-pulse determination. 
The transition between “Type 1”, with a single main 
component having a trailing-side shoulder and a more 
quickly rising leading edge, and “Type 2”, with the main 
component having a leading-side shoulder and a nub fea- 

http://www.cv.nrao.edu/~sransom/presto/ 

We used the lEatough et al.l (I2013blf DM of 1778 cm“^ pc, 
and compared the results to those from times-series dedispersed 
at 0 cm“^ pc and twice the nominal DM in order to discriminate 
between transient RFI and candidate pulses. 


ture on the trailing side, happens more than three weeks 
after the X-ray burst on MJD 56407 and more than two 
weeks before the burst on MJD 56450. The “Type 1” 
shape was seen as ea rly as a week afte r the discovery 
( Eatough et al.l l2013a[) and published in lEatough et al.l 
(|2013b[) . Similarly, the transition between “Type 2” and 
“Type 3”, which has a larger two-peaked component, 
happened more than two weeks after the burst on MJD 
56450 and more than three weeks before the burst on 
MJD 56509. Eor these reasons, we do not associate the 
profile types (which are most likely not absolutely dis¬ 
cretized) with the observed X-ray bursts. Within a single 
observation, the profile shape did not change between 5 
min subintegra ti ons. 

iLvnch et al.l (|2014D also documented the time- 
variability of J1745—2900’s X-band profile as seen with 
the GBT. As their first observation is coincident with our 
last observation, they have also seen the “Type 3” shape, 
which persists and evolves during most of what they have 
labeled a “stable state”. This “stable state” is charac¬ 
terized by relatively smooth profile transitions, a gradual 
flux evolution, and a phase-connected timing solution — 
all in contrast to what they call an “erratic state”, which 
is onset sometime after MJD 56682. Later in their obser¬ 
vations, during epochs with MJDs 56794 and 56865 (both 
in the “erratic-state”), they see a profile resembling what 
we have labeled “Type 1”. We note that we did not wit- 
ness any of t he ver y sporadic profile variability seen in 
iLvnch et al.l (|2Q14D associated with the “erratic state” 
(e.g., the drastic profile changes seen in their last two 
observations, separated by only eleven days), but rather 
we observed each of these three types only for a single 
interval of time. 

3.3. Timing 

Between having bursts, glitches, unstable profiles, and 
timing noise, magnetars are notoriously some of the hard¬ 
est pulsars t o time (cf. the o riginal radio magnetar 
J1809—1943 (iCamilo et al.l 2QQ6h. or see a recent review 
of magnetars in iMereghetti (1201311). As is evident from 
the X-ray and rad io ti ming in iGoti Zelati et all (j2015f ). 
ILvnch et al.l (j2014f ). and iKaspi et al.l (j2014l ). obtaining a 
single phase-connected timing solution for J1745—2900 is 
difficult, due to a significant level of timing noise. Here, 
we measure an overall average spin-down for the purpose 
of summing the data in each epoch. 

Pulse times-of-arrival (TOAs) were measured by cross- 
correlating the time- and frequency-averaged data pro- 
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files with smoothed, “noise-free” template profiles using 
standard PSRCHIVE routines. The templates are gen¬ 
erated by arbitrarily aligning and averaging all of the 
data for which the template is used. Single templates 
were used for the L- and S-band data, but three sepa¬ 
rate templates were used for X-band, depending on the 
profile observed, as discussed in 113.21 Arbitrary phase 
offsets were fit between TO As measured from all of the 
different templates as part of the timing models. These 
phase offsets serve to align the template profiles, but do 
so indiscriminately with respect to pulse broadening from 
interstellar scattering; this has the effect of biasing DM 
estimates if one tries also to measure the dis persiv e delay 
between TO As of different frequencies. See ^3.5.2l for our 
DM measurements based on wideband modeling of the 
L- and S-band data. 

Figure |4] shows the measured values of the spin fre¬ 
quency / as a function of time. The average measured 
spin-down of favg = —8.3(2) x 10“^^ Hz s“^ was suffi¬ 
cient to average the data in each ep och w ith negligible 
smearing for the fiux measurements f 93.4p . and is a rea¬ 
sonable app roximation for the overall trend in the spin 
evolutioirn. This average value also lies between the two 
/ values presented in Table 2 of iKaspi et al.l (j2Q14D for 
the same range of dates. 

Although we are not interested in a full timing solu¬ 
tion for these data in this work, we foun d corrobora- 
tive results when following the suggestion in iKaspi et~al] 
(|2Q14D that there is an abrupt change in / around the 
time of the Sm/i^/BAT-observed X-ray burst on MJD 
56450. Namely, while a single, predictive timing solu¬ 
tion was not found, our pre- and post-burst TOAs are 
described by two simple phase-coherent solutions with 
parameters /pre = —5.005(1) x 10“^^ Hz s“^, /post = 
-9.4799(5) X 10-13 Hz s”!, and /post = -2.696(6) x 
10-^1* Hz s"^. These value s are in good agreement wit h 
thos e in iCoti Zelati et all (|2015f ) . Ikaspi et ^ (|2014D . 
and iRea et al.l (|2013[) . although we were not sensitive 
to /pre- We could only obtain a single phase-connected 
timing solution for all of the TOAs by using five spin fre¬ 
quency derivatives, which is not a predictive ephemeris. 

3.3.1. Profile Alignment 

In Figure O we present the absolute alignment between 
the Chandra 0.3-8 keV X-ray profiles and the GBT radio 
profiles in L-, S-, and X-bands. We determined an inde¬ 
pendent ephemeris for each of the two epochs from the 
radio data by fitting TOAs from each day for only the 
spin frequency, fixing the spin-down parameter at the av¬ 
erage value reported above. These TOAs were measured 
from the frequency-averaged data with 5 min subinte¬ 
gration resolution. The phase-zero time was referenced 
to the arrival of infinite-frequency radiation at the Solar 
System barycenter, which assumes a constant dispersion 
measure of 1778 cm-^ pc between the two observations. 
The X-ray photon arrival times were barycentered also 
using the sky position given in ^2.11 and the JPL Plan¬ 
etary Ephemeris de- 405. These events were folded into 
pulse profiles with 64 phase bins using the correspond¬ 
ing epoch-specific ephemeris by the prepfold program of 

Quantities in parentheses represent the Icr uncertainty on the 
last digit in the respective measurement throughout the paper. 



Figure 4. Average spin evolution of J1745—2900. The three ver¬ 
tical dotted lines correspond to the three X-ray bursts detected 
by Swift/BAT. The two vertical grey bars cover our Chandra ob¬ 
servations. Measurements from the two early S-band observations 
are not included, nor from the X-band epoch on MJD 56425, as 
they were very significant outliers. The quoted uncertainty does 
not include the residual scatter. 

Presto. The alignment based on folding using a single 
ephemeris for both epochs — either the post-burst or the 
multiple frequency derivative ephemeris — yielded indis¬ 
tinguishable results. This is reasonable, since the RMS 
timing residual from either of those ephemerides is on the 
level of individual bins. On the other hand, it may be 
surprising that there seemed to be no interruption in the 
“post-burst” ephemeris; the third detected X-ray burst 
occurred at the midpoint between the two simultaneous 
radio/X-ray epochs. 

We modeled each of the two X-ray profiles with four 
Gaussian components to measure the relative offsets with 
respect to the radio profiles. The offsets and their un¬ 
certainties were determined from Monte Garlo trials, 
where “offset” here refers to the phase that ma ximizes 
a cros s-correlation such as the one prescribed in iTavlorl 
(jl992f ). There was a small offset between the X-ray mod¬ 
els, ^0.02 rot. A difference in DM would shift the rel¬ 
ative phase between the X-ray profile and the S-band 
profile (our fiducial profile) only by ^3x10“^ rot per 
unit DM [cm“3 pc]. Even for the DM difference of 
^17 cm“3 pc measured between these epochs (see 93.5.21 
and Figure [8]), the phase difference is ^0.005 rot. The re¬ 
maining offset can be explained by a combination of the 
variability of the X-ray profile and timing noise, with 
the former being dominant. After removing this differ¬ 
ence, the offsets with respect to the radio profiles do 
not change between the two days within the variance 
of the measurements. The phase offset relative to the 
S-band profile is approximately 0.15(1) rot. The radio 
magnetars J1809—1943 and J1550—5418 both also show 
rough align ment of pulsed radio emission with their X - 
ray profiles (jCamilo et al.l l2QQ7bl : iHalpern et al.l l2QQ8f ) , 
whereas no p ulsed X-ray emission h as been detected from 
J1622-4950 (lAnderson et a,lJl2(m . 

The two double-peaked X-ray profiles appear essen¬ 
tially featureless. The RMS pulsed fractions are given in 
Table [2l There are not sufficient data to decompose the 
profiles into energy bands to look for meaningful spec¬ 
tral dependencies, although we wish to point out a pos¬ 
sible transient feature that appears in the XMM-Newton 
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data recently published bv iCoti Zelati et ahl (|2Q15f ). In 
the e nergy-dependent XMM -Newton profiles of Figure 4 
from ICoti Zelati et alJ (|2Q15[) . there is a conspicuous nar¬ 
row feature on the leading edge of the double-humped 
X-ray profile that is close to the phase of radio emis¬ 
sion (within ^0.05 rot). It appears most prominently 
around phase 0.55 in the 0.3-3.5 keV profile of the third 
XMM-Newton observation (with Obs. ID 0724210501). 
It is also seen in two of the other three energy-dependent 
profiles (except for the highest energy 6.5-10.0 keV pro¬ 
file), contributing to the integrated flux in the energy- 
averaged profile. A similar feature is seen at the same 
phase in the first XMM-Newton observation (with Obs. 
ID 0724210201) to a lesser extent. According to their 
table, these observations were separated by 23 days, 
with the first occurring 19 days after the Chandra ob- 
servations presen t ed he re (which are also included in 
ICoti Zelati et al.l (|2015D b The three Chandra obser¬ 
vations and the one XMM-Newton observation taken 
during these 23 days show no obvious feature, despite 
covering the same range of energies, although Chandra 
recorded only between 10 and 50% of the counts as did 
by XMM-Newton. Therefore, without additional obser¬ 
vations, it remains only a peculiarity. 

3.4. Radio Flux Density 

From the radio data, we made measurements of 
J1745—2900’s flux density as a function of time and fre¬ 
quency. We measured the mean flux densities in 50 MHz 
wide channels and used a weighted average of these mea¬ 
surements to obtain representative flux densities for each 
band, per epoch. 

For all of the L- and S-band profiles, we defined “on- 
pulse” regions as follows. A model pulse profile for each 
frequency was determined from the wideband modeling 
described in SSI We then found the smallest range of 
pulse phases that contained 99% of the integrated flux 
density of the model profile. Examples of the on-pulse 
windows for the scattered L- and S-band profiles can be 
seen in Figure [TJ The mean flux density was calculated 
by averaging the observed flux density in the window and 
scaling it by the duty cycle. The uncertainties were esti¬ 
mated by measuring the mean noise level in the last quar¬ 
ter of each profile’s power spectrun0. We accounted for 
systematics in the residual profile by adding the scaled, 
residual mean flux density to the uncertainty in quadra¬ 
ture. These corrections were small, as the reduced val¬ 
ues of the residuals were usually <1.5 and always <2.0. 

The measurement of the X-band flux densities was 
complicat ed b y the dynamic baseline variations men¬ 
tioned in ^2.11 as well as the intrinsic variability of the 
profile shape. We used polynomial functions to remove 
the baseline variations on a profile-to-profile basis (e.g., 
see Figure [2]). For these profiles, we first centered each 
profile to be near phase 0.5 to avoid edge-effects of the 
polynomial fit from affecting the on-pulse region. A high 
degree polynomial function was fit to the baseline of each 
profile, where in the first iteration an on-pulse window 
with a duty cycle of 6% was blanked out from the fit 
to avoid initially over-estimating the nois4H3- The level 

This is a robust method to estimate the off- pulse variance, 
assuming the profile is resolved (e.g., see IDemorestlDOO’^ i . 

None of the profiles had a smaller duty cycle than 6% and 


of the residual off-pulse noise was calculated, and then 
the on-pulse window was widened until the flux density 
at the edges of the on-pulse region dropped below the 
noise level. The baseline polynomial was then refit to 
the original profile, but with the new on-pulse window 
blanked out. The mean flux density and its uncertainty 
were calculated in these baseline-removed, on-pulse win¬ 
dows as described for the lower frequency data above, 
but a systematic error was added in quadrature to the 
uncertainty that represented the mean flux density across 
the on-pulse phase window removed by the polynomial 
fit. This tested method gives dependable, conservatively 
estimated X-band flux densities. 

3.4.1. Flux Evolution 

The radio flux evolution of J1745—2900 is shown in 
Figure [6l The mean X-band flux density increases 
rapidly in the first half of our observations, increasing 
by at least a factor of ^6 over fifty days, and then ta¬ 
pers off at the 1 mJy level. The earliest reported mea¬ 
surement of J1745—2900’s X-band flux density was ^0.2 
mJy, taken with the Effelsberg 100-m Radio Telescope, 
consistent with our C BT measurement two days later 
(|Eatongh et al.|[2013a[) . Our data show a similar increase 
in the low frequency flux densities. The S-band flux in¬ 
creases by about an order of magnitude over ninety days, 
and in our last five observations covering about thirty 
days, the average L- and S-band fluxes increase by a fac¬ 
tor of two. Given t he mea sured scattering timescales 
for J1745—2900 (see ^3.5.ip and the recently measured 
proper motion of the pulsar, the timescale for refractive 
scintillation to be import ant is much larger th an the span 
of our observations (see iBower et al.1 (|2015f ) for further 
discussion). 

Having picked up where we left off, iLvnch et al.l (j2014l ) 
increased the cadence of CBT X-band observations after 
MJD 56516 and found a similar, slow increase of the 
flux, up to ^3 mJy, over the next 170 days. As al¬ 
ready mentioned, after this “stable state” of slow, steady 
flux increase, the authors found that J1745—2900 en¬ 
tered an “erratic state”, characterized in part by a larger 
and highly variable X-band flux, similar to what was 
seen in two other radio magnetars (jCamilo et al.ll2007bl: 
iLevin et~ani2012[ ) . Superimposed on top of this radio flux 
evolution is a relatively slow decay of the X-ray flux, com¬ 
pared to other magnetap (iRna et al.11201^ : iKaspi et al.l 
I20l4 iLvnch et al.l (20l4 ICoti Zelati et al.l l2015[ ). Be- 
tween our two simultaneous CBT/ Chandra observations 
separated by ^15 days, the radio flux increased by ^60% 
while the X-ray flux decreased by ^10%. This trend 
(seen here and in lLvnch et al.l (j2014f )) is opposite to those 
of the other radio magnetars, which show decreasing ra¬ 
dio a nd X-ray flux wit h time over the course of an out¬ 
burst [Re^^L^ (|2Q12[ ). 

3.4.2. Radio Spectral Index 

Because we have essentially simultaneous observa- 
tion0 of J1745—2900 in frequency bands spaced by two 

a polynomial of degree 15 was used; this was the smallest degree 
polynomial that reasonably and automatically removed systematic 
baseline trends from all of the profiles without having to also vary 
the degree of the polynomial on a profile-to-profile basis. 

In one case, the X-band observation was taken a day earlier; 
see Table [H 
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Figure 5. Absolute phase alignment of J1745—2900’s radio and X-ray profiles determined separately on two days. Note that the brightest 
radio profile is seen in S-band (see Figure [T]). The profiles have 1024 and 64 phase bins, respectively, and are shown as they would be 
observed at the Solar System barycenter for phase-zero MJDs 56499.98000761 and 56515.96999979, refer enced to infinite f r equen cy. The 
assumed dispersion measure is 1778 cm“^ pc. During two later XMM-Newton observations (presented in ICoti Zelati et aTI (1201517 ). there 
is a peculiar, narrow feature seen in the otherwise broad X-ray profile near the phase of radio emission as shown here (see text). 


octaves, we can measure the spectral index 7, where 
Sj^ (X for flux density Si, at frequency v. The up¬ 
per panel of Figure [6] shows 7 as measured between the 
average X-band flux density and the combined average 
flux densities of the lower frequency band(s). The error 
bars were approximated by varying the average fluxes 
within their measurement uncertainties. The decor¬ 
relation bandwidth for diffractive scintillation is much 
smaller than even our native frequency resolution and 
will not be a source of variability here. 


There is no large, obvious stochasticity, as op 

posed 
), but 

to. for examole. J1809—1943 (ILazaridis et al. 

to 

0 

0 

00 

there mav be a trend. IShannon & Johnston 

(120131) re- 


port two early measurements of 7 across the bands span¬ 
ning 4.5-8.5 GHz and 16-20 GHz. The first measure¬ 
ment on MJD 56413 is close to —1.0 in the high fre¬ 
quency band, though it is closer to 0.0 in the lower fre¬ 
quencies, and the second on MJD 56443 is ^ —1.0 across 
both bands, consistent with our measurements more than 
two weeks prior. Our three later measurements indi¬ 
cate a significantly steeper spectrum. The average value 
for 7 of —1.4 is tantamount to the average spectral in¬ 
dex for nor mal pulsars a cross gi gahertz frequencies as 
rep orted in iBates et al.l (20131 ). (Camilo et al.l (j2007c[ ) 
and I Anderson et~al] (2012) both make mention of a gen¬ 
eral steepening of the spectral indices of J1809—1943 
and J1622—4950, respectively, despite remaining much 
fiatter than what is seen in J1745—2900. However, 
(|Lazaridis et al.|[2QQ^ finds the opposite for J1809—1943 
in later observations. 

This finding apparently breaks the mold set by the 
other three radio magneta rs, which have essentially 
fiat (or inverted ) spectra (iCamilo et al.l I2006L 120081: 
iLevin et all 120101: iKeith et ahl 1201ID . However, no firm 
conclusions can be drawn from this handful of measure¬ 
ments from early times in J1745—2900’s outburst, es¬ 
pecially knowing that the other r adio magnetars also 
show a variable radio spectrum (|Gamilo et al.l 12007(4 


iLazaridis et al.l [20081: lAnderson et al.l 120121). In fact, at 
the time of writing, the findings of iTorne et al.l (j2015f ) 
suggest that at much later times (a year after the present 
observations), the radio spectrum of J1745—2900 be¬ 
tween 2 and 200 GHz was much fiatter, with 7 = —0.4(1). 

3.4.3. Spectral Shape 

One example of J1745—2900’s radio spectrum is shown 
in Figure [3 the spectra from the other days are qual¬ 
itatively similar. The spectrum shows a non-power- 
law increase in flux between 1.4 and 2.4 GHz, with a 
possible peak near 2 GHz. The inverted log-parabolic 
shape is reminiscent of what hav e been called “gigahertz- 
peaked spectra” (GPS) pul sars (|Kiiak et al.l 1201 ll 120131: 
iDembska et al.l [2014 I2015F although the GPS pulsars 
supposedly have a much broader spectral shape, over a 
dex in frequency. For reference, we fit a log-parabola to 
the low frequency points, the parameters of which are 
given in the figure. 

It is difficult to explain the spectral shape we see in the 
lower frequencies. It is conceivable that the dense, unique 
environment near J1745—2900 in the Galactic Genter sig¬ 
nificantly alters the spectral shape of radio emission be¬ 
tween 1 and 10 GHz (e.g., via free-free absorption, al¬ 
though the detection of Sgr A* at 330 MHz imp lies a 
low free-free optical depth of <1 (|Nord et al.ll2004f) ). but 
it is difficult to draw any conclusions without a dedicated 
set of observations. 

Another possibility is that we have systematically 
under-estimated the flux: one well known source of bias 
comes from under-estimating the flux at low frequencies 
due to significant area in the scattering tails being lost 
in the calculation of the baseline flux. However, even at 
1.4 GHz the scattering timescale is ^500 ms ~ 0.13 rot 
(see ^3.5.ip . In the worst case of a Kolmogorov scatter¬ 
ing index (—4.4), the scattering timescale at our lowest 
frequency is no more t han 7 ^20% of a rotation. As men¬ 
tioned in iKiiak et~al] (|2011[) and treated graphically in 
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iMacanart et al.l (|2Q1Q[ ) , the pulsed fraction drops by only 
^10% when the scattering timescale is half the pulse pe¬ 
riod. Therefore, we can suggest that at worst we are 
underestimating the L-band flux densities at the ^10% 
level, but this still would imply a positive or approx¬ 
imately flat spectral index between 1.4 and 2.4 GHz; 
the observed flux density drops precipitously somewhere 
thereafter. 

A more promising, albeit provisional possibil¬ 
ity h as been offered up by rece nt modeling of 
the iShannon fc Johnst^ (j2Q13f ) observations. 
iLewandowski et ^ (j2Q15bf ) make a case study of 
J1745—2900 to demonstrate the possibility of thermal 
free-free absorption as the explanation for the GPS. 
For J1745—2900, the authors suggest a combination 
of an expanding ejecta and/or an external absorber 
to explain the cha nging spectrum seen early a fter the 
initial outburst in iShannon fc JohnstonI (j2013f ). The 
free-free absorbed model spectra offer a reasonable 
explanation for the lack of low-frequency detections 
of J1745—2900 immediately after the initial outburst 
and detections above 4 GHz; our two early S-band 
observations may support this idea. Our spectra from 
three months later may also inform the story of an 
evolving or endemic free-free absorbing medium in the 
environment of J1745—2900. 

3.5. Wideband Portrait Model 

As is evident from Figured! J1745—2900 has a highly 
scattered, simple profile across a gigahertz bandwidth, 
from 1.4 to 2.4 GHz. For a nominal DM value of 
1778 cm“^ pc, there is a delay of ^0.66 rotations across 
this band, which is easily measurable. All of the aver¬ 
age L- and S-band profiles showed prominent scattering 
tails from multipath propagation through the interstel¬ 
lar medium (ISM). The quality of the data permitted us 
to make “wideband” measurements of both the DM and 
the scattering timescale r, as well as its power-law index 
a, on an epoch-to-epoch basiJ^. 

For this, w e used the methods an d augmented software 
described in iPennncci et al.l (j2Q14l ) to make a wideband 
“portrait’s model for each of the five epochs where we 
have both L- and S-band observations. For each of these 
epochs we combined the data from the two low frequency 
bands in a fit for a global portrait model that included a 
single scattered Gaussian component with profile evolu¬ 
tion parameters, a constant baseline term, a phase offset 
between the bands, DMs for each band, and the scat¬ 
tering index a. The scattering timescale is defined in 
the usual way by assuming a one-sided exponential pulse 
broadening function for the ISM, so that an observed 
profile p{(p) is the convolution given by 

p{ip)= ( 1 ) 

where p is the rotational phase, is the spin period, 
H is the Heaviside step function and g{ip) is the intrinsic 

The two earliest S-band observations were exceptions; cor¬ 
rupted data, low signal-to-noise ratios, and the lack of L-band 
data resulted in uninformative measurements of the DM, r, and 
a. This is also the reason thes e ob servations were excluded from 
the average / measurement in iS31 

We use the word “portrait” to mean the total intensity profile 
as a function of frequency. 


total intensity profile shape. For a power-law spectrum of 
density inhomogeneities in the ionized ISM r is expected 
to have a power-law dependence on frequency n as 

= (2) 

with reference frequency In all cases, the scatter¬ 
ing timescale (133 ms at 2 GHz; see below) dominates 
the smearing from the process of incoherent dedispersion 
(^0.7 ms (z//2 GHz)“^), the smearing from an incorrect 
DM when averaging channels ( r\j 25 yus ((5DM/cm ^ pc) 
(z//2 GHz)“^), and the temporal resolution (1.8 ms for 
2048 profile bins), so we have not included those modifi¬ 
cations of the pulse profile shape in the model. However, 
deviation from the simple timing models discussed in ^3.31 
(e.g., see Figure |4]) during any of these epochs could add 
profile smearing in the integrated profiles (at the level 
of ^tens of ms — a significant fraction of the scattering 
timescale). We avoided this source of bias by iterating 
over the timing model to remove the timing residual on 
a per-epoch basis. 

We model g with a frequency-dependent Gaussian 
function, 

g{v, (fi) = A(y) exp ( - 41n(2) 

which is parameterized by its location full-width- 

at-half-maximum (FWHM) (7(z^), a nd am plitude A{n). 

As described in IPennncci et aid (|2Q14[) . each of these 
parameters nominally has an additional parameter de¬ 
scribing its frequency dependence. However, because 
this combined band has a fractional bandwidth of “only” 
^0.5, we assume (Pg{n) = (/?o is a frequency-independent 
value. That is, we assume there is no drift intrinsic to 
the one component across the band. 

Furthermore, when allowing for a frequency-dependent 
cr, we found no significant evolution, and so we chose 
also to fix the evolution cr(z/) = do to be a frequency- 
independent fit parameter in our final portrait models. 
This choice was further justified by performing indepen¬ 
dent per-channel profile fits of a single, scattered Gaus¬ 
sian component and examining the frequency evolution 
of a. Also, there are X-band observations for three of 
these epochs, and in these cases the FWHM of the X- 
band profiles (all of “Type 3”), when fitted with a single, 
unscattered Gaussian component, was always within the 
scatter of those measured from the lower frequency ob¬ 
servations. These results are consistent with the weak (or 
lack o f) frequency dependence of a found in iSpitler et al.l 

(l20a . 

We normalized the intensities of the data to be fit by 
the maximum profile value in each frequency channel to 
remove the unusual spectral shape (see ^3.4.3p . This al¬ 
lowed the Gaussian amplitude to be easily modeled by a 
power-law function for A{u). In all cases, the reduced 
of the fit was <1.1, and a second Gaussian component 
was never justified by the residuals. 

The combination of the quality of the X-band data, the 
variability of the profile, and the expected value of r at 
8.9 GHz (<1 ms, comparable with our native time reso¬ 
lution) was such that we did not attempt to incorporate 
this high frequency data into our wideband profile model, 
nor did we measure the scattering timescale in either the 
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Figure 6. The early radio flux (bottom panel) and spectral (top panel) evolution of J1745—2900 over 100 days from the observations 
in Tabled The vertical demarcations are the same as in Figure 0] I Lynch et al.l (|2014l ) And a continuation of the slow, steady increase 
in X-band flux for another six months, which is followed by what they call an “erratic state”. The apparent excess average S-band flux 
density during MJD 56501 is explained by the fact that the lower half of the band was corrupted (see Table [T}, and the pulsar’s flux density 
apparently increases with frequency in this range (see Figure [7|. The average value of the spectral index 7 is about —1.4; see text for 
details. 



Figure 7. An example of J1745—2900’s radio spectrum from the 
brightest observed epoch, MJD 56516. The markers are as in Fig¬ 
ure [ 6 ] note that the flux densities agree in the ~100 MHz overlap 
between L- and S- band. A similar inverted parabolic shape over 
log-frequency is seen during the other sets of (nearly) simultane¬ 
ous observations, which is reminiscent of the so-called GPS pulsars. 
The coefficients a, 5, and c of the fitted dashed parabola (logio(Si^) 
= ax^ + 6 a: + c, for x = logio(^')) are given in the plot, along with 
the spectral index 7 , which for this plot was fitted between the 
peak of the parabola and the X-band data. 

av erag;e profile o r the singl e pulses. We r efer t he reader 
to iBower et al. ] (|2Q14f ) and iSoitler et al.l (|2Q14D for high 
frequency scattering measurements of J1745—2900. 

3.5.1. Pulse Width & Scattering Parameters 

The results from our wideband models are shown 
in Figure [8l There was no significant change 
in the measured FWHM of the unscattered pro¬ 
file, and our average (frequency-independent) value of 
91.9(4) ms = 0. 0244(1) rot is also c onsistent with what 
is reported in iSpitler et al.l (j2014f ). The scattering 
timescale at 2 GHz, r2GHz, appears to increase by ^10% 
over the four weeks, and the scattering index a deviates 
from its average value, first to a Kolmogorov value near 
—4.4 (the dash-dotted line), and then to a much shal¬ 
lower value near —3.0. Both of these results are some- 
w hat peculiar, but sim ilar variations are also reported 
in ISpitler et al.l (|2014D , though they do not discuss the 
temporal evolution of either quantity. That is, their pub¬ 


lished values of r from a variety of epochs and frequencies 
cannot be unified by a single scattering timescale and in¬ 
dex. In fact, their measurements of r show more scatter 
over the course of their observations than those presented 
here, which have some overlap. When the authors com¬ 
bine all of their measurements, they find an average value 
for G of -3.8(2). 

We checked our measurements in two ways. First, we 
performed conventional profile fits of a single, scattered 
Gaussian component to each individual frequency chan¬ 
nel, independent of any evolutionary constraint. The 
values of a, r, and a for each epoch were consistent 
with what we found by applying the wideband model¬ 
ing method. In Figure [9l we show the measurements of 
r measured in this way for the brightest observed epoch 
(MJD 56516) and over-plot the fitted power-law, which 
has the most extreme a value of the five epochs. There 
was nothing unusual about the data from this epoch 
in terms of RFI, data removal, calibration, or baseline 
variations. Second, as a check for our average values, 
we summed all of the data portraits together by coher¬ 
ently stacking the observations (having fit for a phase 
and DM in each epoch), and fit a single wideband model 
to the averaged data (with the same constraints as ear¬ 
lier). Using this method, we obtained similar average 
values: a = 0.0246(1) rot, r2GHz = 133.0(5) ms, and 
a = —3.71(2), the la tter of which i s in c oncert with the 
average a value from ISpitler et al.l (j2014f ) . Our extrapo¬ 
lated value of tighz = 1.74(3) s is only slightly at odds 
with their average value of tighz = 1.3(2) s, which is 
probably due to th e temporal variabil ity of r. As oth¬ 
ers have noted (e.g.. lBower et al.l (j2014l )). the anticipated 
value for tighz along this line of sight based on empirical 
relationships, for a DM of 17 78 cm~^ pc, is about 600 x 
larger than what is observ ed (jKrishnaknmar et al.ll20l5 
iLewandowski et al.ll2015a[ ). 

iBower et al ] (l2nTl determined ts.tghz ^ 2 ms from 
interferometric measurements of J1745—2900’s single 
pulses, implying that a value for a as shallow as —3 
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Figure 8. Our wideband measurements from five epochs. The 
vertical demarcations are the same as in Figure [4] The FWHM 
showed neither frequency nor temporal dependence. The trend 
in the scattering timescale r is less scattered and more precise 
than the measurements presented in ISoitler et al.l (|2014l b The 
dash-dotted line in the panel for the scattering index a marks the 
fiducial Kolmogoro v value of —4.4, and the dashed lines mark the 
ISoitler et al.l (120141 ) measurement of —3.8(2). The additional mark¬ 
ers in the bottom panel are the same as in Figure[6l the blue/down¬ 
pointing triangles are L-band measurements, and the green/right¬ 
pointing triangles are S-band measurements — the dots are their 
weighte d average. The dashed lines here are the lEatough et al.l 
(I2013bl) DM of 1778(3) cm“^ pc. See the text for a discussion of 
the DM measurements. 


is not unbelievable. Furthermore, scattering measure¬ 
ments from two high-DM pulsars discovered near the 
Galactic Center (both within 0.3° and having DMs 1100- 
1200 cm“^ pc) implied a = —3.0(3) (|Johnston et al.l 
Hooi. It is not uncommon for pulsars to have a > — 4, 
particularly along special lines-of-sight, and it is em¬ 
pirically suggested that the highest DM pulsars may 
have an a verage scattering index significantly shallower 
than — 4 (jLohmer et al.l 1200 iL l2004 iLewandowski et ^ 
l2015a|) . Note that observing a ^ —4.4 does not nec¬ 
essarily imply a non-Kolmogorov spectrum of density 
inhomogeneities; rather, it could be t hat a non-thin- 
screen geometry may be resp onsible (|Cordes fc Laziol 
I 2 OOII: ILewandowski et ^120131 ). 

3.5.2. Dispersion Measures 


of the Galactic Center Magnetar 


11 


10 ^ 


Figure 9. Independent per-channel measurements of the scatter¬ 
ing timescale r and the fitted scattering index a for the brightest 
set of L- and S-band observations, on MJD 56516; similar plots 
from the other four days have a significantly more negative slope. 
The dashed lines represent our measurement, whereas the dotted 
lines show our average value of a from wideband modeling and 
the fiducial Kolmogorov value of —4.4, all with the same value of 
r2GHz- Here, the measurement uncertainties have been infiated by 
the reduced 2. 

The bottom panel in Figure [8] shows the best-fit 
DMs as determined by the wideband models in the es¬ 
sentially simultaneous L-band and S-band observations 
(blue/down-pointing and green/right-pointing triangles, 
respectively). The black points are the weighted average 
of the two measurements; there is obviously some vari¬ 
ance about the nominal value of 1778(3) cm“^ pc, and 
our overall average value is ^1781(1) cm“^ pc. Without 
exception, the measured L-band DMs are greater than 
those measured in S-band. There is also only one epoch 
where the la uncertainties have any overlap; the RMS 
variance of the differences is ^6 cm“^ pc. The absolute 
DM differences cause residual dispersion on the order of 
<20 ms ^ 10 bins (for 2048-bin profiles) across the cor¬ 
responding band, and so they present significant profile 
deviations. To make sense of the discrepant DMs be¬ 
tween the two frequency bands, we consider a number of 
possibilities. 

First, the time-averaged data for each epoch showed 
few systematics with negligible baseline variations, so we 
do not believe that data quality was an issue here. 

Next, as is well known, the measured absolute DM 
will be affected by the choice of profile alignmenlF^. We 
can rule out any simple, constant profile evolution as 
the source of the differing DMs because such a mod¬ 
ification introduces a constant difference in the DMs; 
the changes in the measured DM should be the same 
independent of the choice of alignment. Even if our as¬ 
sumption that there is no intrinsic drift in the location 
of the (unscattered) profile component across the band is 
wrong, allowing for a drifting component will still repro¬ 
duce discrepant DMs; we have verified this by allowing 
for frequency evolution in the location parameter of the 
Gaussian component. 

A second confounding element from our modeling 
could be the use of different models for each epoch; if 
they are all systematically wrong in their alignments or 

For example, DMs are significantly biased when either as¬ 
suming a constant profile shape in the presence of scattering, or 
aligning scattered profiles by conventional methods because the 
convolution of the ISM pulse broadening function with a profile of 
finite width introduces a delay that is a function of the scattering 
timescale (i.e. frequency). This is partly why proper wideband 
modeling is necessary. 
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representation, they could be wrong differently. One way 
to check this is to simply use one fixed model to remake 
the DM measurements. We used the average portrait 
model discussed earlier and confirmed that the DMs re¬ 
main similarly extreme, within ^2 cm“^ pc, comparable 
with the measurement uncertainties. In fact, we tried 
a large number of fixed and variable portrait models, 
but never obtained either consistent DMs or DMs with 
a near constant offset. So, to the extent that r and/or 
a are measurably changing, we are justified in keeping 
them as free parameters for each epoch’s model. 

Similarly, the known profile variability that is seen in 
all of the radio magnetars could also play a role when 
using either a fixed or variable portrait model. However, 
besides the fiux density, any underlying profile shape 
changes either with time or frequency are masked by 
the large level of scattering. As mentioned, the FWHM 
does not seem to change significantly in either time or 
frequency. Furthermore, the three X-band observations 
taken during these epochs show no large profile changes, 
and are all of the “Type 3” shape. 

Next we can ask whether or not the slight asynchronic¬ 
ity could have any effect; that is, could the DM change 
so significantly on ^hour timescales? We will return to 
this question below, but it is not an uncommon a priori 
assumption to expect that the observed DM does not 
change between observations separated by < 4 hr. 

One could also ask if the method by which we mea¬ 
sure the DMs introduces a systematic error, where the 
error may depend on the exact values of r and a, or 
even the spectral shape. To answer this, we performed 
a number of Monte Carlo simulations. In the simula¬ 
tions, we used the models from MJDs 56500 and 56516, 
which have the most extreme values for the difference 
between the DMs, and the most extreme values for r 
and Of, respectively. For each trial, we made fake L- and 
S-band observations by appropriately constructing the 
model for that band and scaling each frequency chan¬ 
nel’s amplitude to match the spectral shape. We then 
added random frequency-dependent white noise to the 
model at the same level as measured from the data por¬ 
traits and finally dispersed the fake data with a DM of 
1778 cm“^ pc. Visual inspection verified that the fake 
data were faithfully rendered. We used the same method 
to meas ure the DM (an d phas e), which is described in 
detail in iPennncci et al.l (j2014f ) . In summary, the mea¬ 
sured DMs were always in accord and unbiased, and the 
uncertainties were accurately estimated. We conclude 
that the measurement method produces accurate DMs, 
independent of the model parameters, provided that the 
model for the data is accurate. 

We assume in our measurements that the phase offsets 
(A0) incurred by finite-frequency signals due to propaga¬ 
tion through the ionized interstellar medium scale as pre¬ 
dicted by the usual cold-plasma dispersion law such that 
A(f) (X This is certainly the case to first-order 

even over large, low frequency bandwidths (jHassall et al.l 
|2Q12[ ). However, to the extent that we understand the 
ISM to be mhomogeneous — after all, we do observe 
pulse broadening — then it is anticipated that the sim¬ 
ple dependence will be an insufficient description at 
some level for broadband DM measurements. When an 
inhomogeneous medium causes multi-path propagation 


of radio waves where the path depends on frequency, the 
sampled column density of free electrons (the DM) will 
also be a function of frequency. Thus, we are left with the 
intriguing possible explanation that the DM inconsisten¬ 
cies we are seeing are the consequence of imposing a 
dispersion law onto a frequency-dependent DM (DM(z/)) 
due to an inhomogeneous ISJVP^. 

To our knowledge, the most recent claim for hav- 
in g observed fr e qnerm y-dependent DMs was reported 
in lAhuja et al ] (l 2 nnl for the slow, low DM pulsars 
B0329-F54 and B1642—03, although they observed lower 
DMs at lower frequencies. However, the authors only 
made one set of simultaneous pairs of dual-frequency 
measurements per pulsar. We argue that to confidently 
segregate the effects of profile evolution, DM variations 
with time (DM(t)), DM(z/), and other potential con¬ 
founding factors, many epochs of simultaneous, wide¬ 
band (large fractional bandwidth) observations of a sta¬ 
ble, preferably high DM pulsar need to be made. A sim¬ 
ilar re commendation was recently made by iCordes et al.l 
([20T^ in their detailed study of frequency-dependent 
DMs, which makes theoretical predictions for the char¬ 
acteristic timescales and sizes of D M(z/) effects. 

In their treatment of the problem, ICordes et al.l (j2Q15l ) 
predict the minimum scale of DM variations about a 
mean value, 

DMrms ^ 0F/Are ^ 3.84 X 10“^ cm“^ pc z/qhz (4) 

where z^ghz is the frequency in GHz and is the size 
of the phase perturbations over the Fresnel scale. If = 
y^[cD)/(27n7)^ for the speed of light c and source distance 
D. For J1745—2900, which is in the strong scattering 
regime, (pp will be very large. We estimate it from their 
prescription, 

~ 9.6 rad (^/^) ^ , (5) 

where Aud is the scintillation bandwidth, which is 
readily estimated from our scattering measurements 
as ^ 1.16/(27rr(z^)). For 1.4 and 2.4 GHz, we find 
DMrms ^ 10 and 5 cm“^ pc, respectively. These can 
be compared to the RMS DM values as measured in L- 
and S-band of ^ 9 and 7 cm“^ pc, respectively. The char¬ 
acteristic spatial size for the DM differences near 2 GHz 
will be several Fresnel scales, which can be converted 
to a characteristic time by usi ng the recently mea sured 
proper motion of 236 km s“^ (| Bower et al.ll201^ . For 
our range of frequencies, the characteristic timescale as¬ 
sociated with the Fresnel scale size is ^3 hr, comparable 
to the separation between the observations on a given 
epoch. Therefore, it may be that the small temporal gap 
between the observations contributes somewhat to the 
difference in the DMs, but we certainly do expect that 
the DMs vary significantly on different days, separated 
by many Fresnel timescales. 

Finally, ICordes et al.l (j2Q15l ) make a prediction for the 
observed RMS difference between DMs at frequencies u 

This is opposed to other supposed origins of DM(z/) relating to 
magnetospheric propagation effects or magnetic sweepback, which 
wo uld likely have differ ent statistics from an I SM induced DM(i/); 
see lHassall et al.l (|2012l ) or lAhnia et al.l (|2007l i for an overview. 
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and z/', 

~ 4.42 X 10“® cm-3 pc F^ir) (y^)> (6 ) 

where we have ignored a geometric factor of order unity 
and the function contains the frequency dependence 
for r = yjy' given the power-law index (3 for the 
wavenumber spectrum of density inhomogeneities. For 
y = 2.4 GHz and = 1.4 GHz, ^ 4 cm“^ pc, com¬ 
pared to our observed RMS difference of ^6 cm“^ pc. 

That the predicted and observed values are similar 
may be coincidence, but we note the corroborating facts 
that J1745—2900 is the highest DM pulsar, is relatively 
bright, highly scattered, has a simple, easily modeled 
profile, and does not show significant profile evolution 
or stochastic profile variability (at least in these obser¬ 
vations). Furthermore, we verified that our measure¬ 
ment method produces inconsistent (and biased) DMs 
between the bands by introducing non-z^“^ phase delays 
into our fake data simulations described earlier. After 
ruling out the other potential sources for the inconsistent 
DMs, we suggest that J1745—2900 may have an observ¬ 
able frequency-dependent dispersion measure. 

A potential counter argument is that over many Fres¬ 
nel timescales, one expects the sign of the DM differences 
to change, such that the observed low frequency DM be¬ 
comes smaller than the high frequency DM. Between the 
small number and low density of epochs, the potentially 
incorrect portrait model, and ISM uncertainties (the pre¬ 
dictions here are based on a thin-screen model with a 
Kolmogorov spectrum of density p erturbations, which is 
partly supported by the findings in iBower et al.l (j2Q14f )). 
it is conceivable that this observation is not inconsistent 
with a frequency-dependent DM as described. 

Determining whether or not a difference in DM as 
seen in two frequency bands is intrinsically a DM(z/) ef¬ 
fect is complicated by the issues described above, and 
with only five measurements we obviously cannot draw 
any definite or statistical conclusions, but future studies 
could potentially disentangle the evolution of DM(t,z/), 
th e profile, and other ISM parameters. One strategy, 
as iCordes et al.l (|2Q15f ) note, is to model the frequency 
dependence of the dispersive delays as something other 
than y~‘^. This should be done for many epochs, at least 
as long as the timescale for refractive scintillation, over 
which time the specific frequency dependence of the aver¬ 
age DM remains stable. For J1745—2900, this timescale 
is potentially many years. 

4. SUMMARY & DISGUSSION 

In this paper we have presented multi-epoch, multi¬ 
frequency wideband GBT observations of the Galac¬ 
tic Genter radio magnetar J1745—2900 at 1.5, 2.0, and 
8.9 GHz from the first ^100 days after it was discov¬ 
ered. After its initial X-ray burst on 25 April 2013, 
J1745—2900 underwent two additional bursts in the 
course of our observations. For two epochs, during which 
time we collected data from three radio bands, we also 
have simultaneous X-ray observations taken with Chan¬ 
dra. An analysis of the radio data, as well as a joint 
analysis with the X-ray data, yielded a few noteworthy 
results. 

1. We found no anomalous radio bursts or giant-pulse- 


like individual pulses in any of the X-band ob¬ 
servations. Similarly, the smooth transitioning of 
the X-band profile between three broad categories 
seems to have also been unperturbed by the Swift- 
detected X-ray bursts. 

2. Our simpl e radio timing anal ysis corroborates the 
findings o flKaspi et al.l (12014), wh ich are also sup¬ 
ported bv iCoti Zelati et al.l (j2Q15f ). We presented 
the absolute alignment of the three radio and 0.3- 
8 keV profiles. The near-alignment of the radio 
components with the X-ray profile is similar to 
the two other radio magnetars that have published 
alignments. We also mak e note of a possibl e tran - 
sient X-ray feature from iCoti Zelati et al.l (j2Q15f ) 
because of its proximity to the phase of radio emis¬ 
sion located ^ 0.2 in phase preceding the peak in 
the X-ray profile. 

3. The evolution of our early radio flux measurements, 
showing a relatively stable growth from around 
the time of the initial outburst, is consistent with 
th e continued GBT X- band observations presented 
in iLvnch et al.l (j2014f ) and with what they have 
called a “stable state” B. The combination of the 
gradual flux evolution with the simple timing and 
profile variability results leads us to extrapolate 
J1745—2900’s “stable state” back to the time of 
its initial burst. 

4. The shape of J1745—2900’s low frequency radio 
spectrum is potentially positive or fiat, whereas 
it shows a “typical” spectral index of ^ —1.4 be¬ 
tween ^2 and 9 GHz, at least during a brief pe¬ 
riod ^100 days after its initial outburst, around 
the times of two later X-ray bursts. This steep 
spectral index might indicate a different mag- 
netospheric configuration during these times, al¬ 
though the evolving spectra may be a result 
of environmental factors a nd free-free absorption 
(jLewandowski et ^l2Q15b[ ). The possible variabil¬ 
ity of 7 means that dedicated observations cov¬ 
ering several higher frequency bands need to be 
carried out over many epochs to confirm this (cf. 
Irbrne et al]l2Q15f ). 

5. We made wideband models of J1745—2900’s low 
frequency radio “portrait” to measure the scatter¬ 
ing timescale, scattering index, and the DM as 
a function of time. Our average measurements 
are consisten t with what has been published in 
iSpitler et al.l (j2014f ). though the ISM parameters 
may be variable. Time-variable scattering param¬ 
eters would complicate the predicted sensitivities 
of future pulsar surveys of the Galactic Genter. 
Lastly, we make a suggestion that our discrepant, 
nearly simultaneously determined DMs are a man¬ 
ifestation of an ISM-induced frequency-dependent 
dispersion measure, and that future observations 

While our observations were taken over a shorter range of 
time (about a third), our cadence of observations is comparable 
to theirs taken during the “erratic state”, the onset of which was 
apparently unrelated to X-ray bursts. 
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could make a case study out of J1745—2900 to in¬ 
vestigate DM(z/) — provided the pulsar remains 
visible and stable. 

J1745—2900 shares several (but not all) properties 
with the other three radio mag netars, J1809—1943 , 
J1550-5418, and J16p -4950 (iCaTnilo et, a], I 
I2nn7al: iTjevin et, alJ 1201(11) . Common properties of the 
pulsed radio emission from magnetars are: a) a delay 
in the appearance of the radio emission after the X- 
ray outburst onset, b) variable pulse profiles and radio 
flux on timescales from hours to days, c) a large rota¬ 
tional (spin-down) luminosity with respect to the qui¬ 
escent X-ray luminosity, d) a decrease of the radio flux 
as the X-ray flux decays, and e) a flat radio spectrum 
over a wide range of frequencies. J1745—2900 grossly 
shares the first three properties with the rest of its class. 
However, while in all other cases the radio flux was ob¬ 
served to decay as the X-ray outburst was fading, the 
long-term radio and X-ray flux evolution of J1745—2900 
is at variance with this trend. The radio flux shows a 
re-brightening hundreds of days after the outburst on¬ 
set and the X-ray emission is decaying very slowly, chal¬ 
l engin g current crustal cooling models iCoti Zelati et al.l 
(Hoill). Furth ermore, the r ecent ly published flux mea¬ 
surements by iTorne et~ld ] (IMsI) taken one year after 
those presented here suggest that the 8.35 GHz flux re¬ 
mained stable at the ^3 mJy level over thirty days. An¬ 
other interesting peculiarity of this radio magnetar was 
the steep (and possibly free-free absorbed) radio spec¬ 
trum seen in o ur observati o ns, th ough the more recent 
observations in iTorne et al.l (j2Q15f ) suggest that the spec¬ 
trum has since flattened. 

Of particular interest is J1745—2900’s low quies¬ 
cent luminosity co mpared to its hig h rotational power 
(Ax,qui/Arot < 1; iRea et aP l2013f ). This peculiar¬ 
ity of the four radio magnetars, which is at vari¬ 
ance with canonical magnetars (for which the fact that 
Ax,qui/Arot > 1 A^s always been u sed as proof of their 
magnetically dominated emission (iMereghetti fc Stellal 

119951: [Thompson fc Duncanlll995l:[Mereghettill2008f )). h^ 

been viewed as evidence for a similar mechanism power¬ 
ing the radio emission from magnetars and normal pul¬ 
sars alike. In fact, while normal radio pulsars have pri¬ 
marily dipolar-dominated magnetic fields (Hp), magne¬ 
tars have a substantial toroidal component (H^) that is 
present in both the internal and external fields. This 
toroidal component is the main reason for their quies¬ 
cent X-ray luminosities, hot surfa ce temperatures, flaring 
emission, and outb urst activity (jThompson et al.l I2QQ2I: 
iBeloborodovI l2QQ^ . For a fixed dipolar field, the inter¬ 
nal toroidal held has no significant effect on the lumi¬ 
nosity unless Bfh > _ , a s is the case for most mag¬ 

netars (jVigano et al.ll2Q13f ). Both radio magnetars and 
high-H radio pulsars have systematically lower toroidal 
fields and higher rotational energies than typical magne¬ 
tars; this is in agreement with the former being fainter in 
quiescence and having a softer X-ray spectrum (a lower 
crustal toroidal held results in less heating prod uced by 
Joule dissipation in the crust, iPons et ffll2QQ9[ ). As for 
the energy powering the radio emission, simulations of 
high dipolar held pulsars that have a small toroidal com¬ 
ponent showed that the particle acceleration and sub¬ 
sequent ignition of the cascade process could proceed 


as it does in normal pulsars, successfully reaching the 
open- fleld line region and generating pulsed radio emis¬ 
sion (jMedin fc Fail IThqI ) . On the other hand, for an 
extremely strong toroidal component, it is expected that 
the particle cascades cannot reach the open-fleld lines 
due to the powerful currents formed as a consequence of 
the twisted magnetosphere. Radio magnetars might lie 
in between, having a high enough rotational energy to 
power pair cascades as in normal pulsars, but also hav¬ 
ing toroidal components lower than typical magnetars, 
resulting in lower quiescent X-ray luminosities. 

In the above picture, the possible radio flux increase, 
the steep spectrum, and the slow cooling of the X-ray 
outburst might be explained by the presence of a strongly 
twisted bundle, which can account for the radio emission 
and the additional heating by particles slamming onto 
the surface. If the radio emission is generated by ac¬ 
celeration of particles only in this part of the magneto¬ 
sphere, then the radio flux and the X-ray flux might be 
unrelated. In particular, untwisting of the bundle during 
the outburst decay might induce fewer currents blocking 
the pair cascade generation, hence more radio emission 
from this region. However, these are only speculative, 
plausible hypotheses. Proof of this scenario would need 
a longer monitoring of the radio and X-ray emission, as 
well as detailed magnetohydrodynamical simulations of 
particle acceleration and pair cascades in a strongly mag¬ 
netized and twisted bundle. 
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